NUMERICAL SCHEMES FOR G EXPECTATIONS 



YAN DOLINSKY 
DEPARTMENT OF MATHEMATICS 
ETH, ZURICH 
SWITZERLAND 

Abstract. We consider a discrete time analog of G-expectations and wc prove 
that in the case where the time step goes to the corresponding values converge 
to the original G-expectation. Furthermore we provide error estimates for the 
convergence rate. This paper is continuation of 0]. Our main tool is a strong 
approximation theorem which we derive for general discrete time martingales. 



1. Introduction 

In this paper we study numerical schemes for G-expectations, which were intro- 
duced recently by Peng (see [7] and [8]). A G-expectation is a sublinear function 
which maps random variables on the canonical space := G([0, T]; R'*) to the real 
numbers. The motivation to study G-expectations comes from mathematical fi- 
nance, in particular from risk measures (see [6] and [9]) and pricing under volatility 
uncertainty (see [5] ,[5] and [H]). 

Our starting point is the dual view on G-expectation via volatility uncertainty 
(see ^), which yields the representation ^ — )■ suppgp Ep[£] where V is the set of 
probabilities on G([0, T]; M'^) such that under any P gV, the canonical process B 
is a martingale with volatility d{B) /dt taking values in a compact convex subset 
D C S^j. of positive definite matrices. Thus the set D can be understood as the 
domain of (Knightian) volatility uncertainty and the functional above represents 
the European option (with reward ^) super-hedging price. For details see ([2] and 

i)- 

In the current work we assume that ^ is of the form F{B, (B)) where F is a path- 
dependent functional which satisfies some regularity conditions. In particular, ^ can 
represent an award of path dependent European contingent claim. In this case the 
reward is a functional of the stock price which is equal to the Boolean exponential 
of the canonical process, and so quadratic variation appears naturally. 

In [4], the authors introduced a volatility uncertainty in discrete time and an 
analog of the Peng G-expectation. They proved that the discrete time values con- 
verge to the continuous time G-expectation. The main tools that were used there 
are the weak convergence machinery together with a randomization technique. The 
main disadvantage of the weak convergence approach is that it can not provide er- 
ror estimates. In order to obtain error estimates we should consider all the market 
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models on the same probability space, and so methods based on strong approxima- 
tion theorems come into picture. In this paper we consider a bit different (than in 
[1]) discrete time analog of G-expectation and prove that in a case where the time 
step goes to the corresponding values converge to the original G~expectation. 
Furthermore, by deriving a strong invariance principle for general discrete time 
martingales, we are able to provide error estimates for the convergence rate. 

The paper is organized as following. In the next section we introduce the setup 
and formulate the main results. In Section 3 we present the main machinery which 
we are use, namely we obtain a strong approximation theorem for general martin- 
gales. In Section 4 we derive auxiliary lemmas that we use for the proof of the main 
results. In Section 5 we complete the proof of Theorems 12.21 - 12731 

2. Preliminaries and main results 

We fix the dimension rf e N and denote by || • || the sup Euclidean norm on W^. 
Moreover, we denote by S'' the space of d x d symmetric matrices and by S!^ its 
subset of nonnegative definite matrices. Consider the space S*^ with the operator 
norm ||y4.|| — sup||^||^]^ ||A(t;)||. We fix a nonempty, convex and compact set D C 
§'j_; the elements of D will be the possible values of our volatility process. Denote 
hy n = C{[0,T];R'^) and T = C{[0,T];S'^), the spaces of continuous functions 
with values in R'^ and S"^, respectively. We consider these spaces with the sup norm 
||x|| = supQ<f<;y ||a;t||- Let F:f2xr^-Mbea function which satisfies the following 
assumption. There exits constants Hi, H2 > such that 

(2.1) \F{ui,vi)-F{u2,V2)\ < ffiexp(i/2(||ui|| + 1^211 + ||fl|| + ||«2||)) X 

(||ui - U2II + - W2II), Ui,li2ef2, Vi,V2eT. 

Without loss of generality we assume that the maturity date T — 1. We denote by 
B = {Bt)o<t<i the canonical process (on the space ft) Bt{uj) — uJt, u d ft and by 
J^t ■— (^{Bs, < s < t) the canonical filtration. A probability measure P on fl is 
called a martingale law if i? is a P- martingale (with respect to the filtration 
and Bq = P-a.s. (all our martingales start at the origin). We set 

(2.2) Vd = {P martingale law on D, : d{B) /dt e D, P x dt a.s.}, 

observe that under any measure P € Vd the stochastic processes B and (B), are 
random elements in ft and F, respectively. Consider the G-expectation 

(2.3) V= sup EpF{B,{B)) 

PeVr, 

where Ep denotes the expectation with respect to P. A measure P G Vn will be 
called e-optimal if 

(2.4) V < e + EpF{B,{B)). 

Our goal is to find discrete time approximations for V. The advantage of discrete 
time approximations is that the corresponding values can be calculated by dynam- 
ical programming. Furthermore, we will apply these approximations in order to 
find e-optimal measures in the continuous time setting. 

Remark 2.1. Let S ~ {{S} , 8^)}]^^ be the Doolean's exponential £{B) of the 
canonical process B, namely SI = S'^exp (BI - {B')t) , i < d, t G [0,1]. The 
stochastic process S represents the stock prices in a financial model with volatility 
uncertainty. Clearly any random variable of the from g{S) where g : G([0, T]; M'') — >■ 
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R+ is a Lipschitz continuous function, can be written in the form g{S) = F{B, (B)) 
for a suitable F which satisfies i2. Thus we see that our setup includes payoffs 
which correspond to path dependent European options. 

Next, we formulate the main approximation results. Let be a distribution on 
R'^ which satisfies the following 

(2.5) / xdv{x) — and / x^x-'dv{x) — Sij, 1 < i, j < d 

where Sij is the Kronecker-Delta. Furthermore, we assume that the moment gen- 
erating function ipi^iv) ■— XcgRd exp(^^^^ a:*2/*)(ii'(a;) < oo exists for any y g R'^ 
and, for any compact set K C R"^ 

(2.6) sup sup ^ ) <oo. 



neN yeK 



In 



Observe that the standard d-dimensional normal distribution v = N{0,I) is satis- 
fying the assumptions (I2.5p - (l2.6p . 

Let n € N and Yi, F„ be a sequence of i.i.d. random vectors with C{Yi) = v, 
i.e., the distribution of the random vectors is i^. We denote by A'^ the set of all 
d-dimensional stochastic process M = {Mq, ...,Mn) of the form, Mq = and 

* 1 

(2.7) M, = ^— (/.,(Yi,...,y,_i)y„ l<i<n 

j=i 

where 0j : (R'')-'^^ \/D := {^/A : A e D} and li, F„ are column vectors. As 
usual for a matrix A G S'J^ we denote by ^/A the unique square root in . Observe 
that Af is a martingale under the filtration which is generated by Yi,...,y„. Let 
(M) be the (§5j_ valued) predictable variation of M. In view of (|2.5p we get 



1 



k 



(2.8) (M)fe = -> ^j(Yi,...,Y,_i), l<k< 



n ^ ■' 



n 



and we set (M)o = 0. Let W„ : (R'')"+i x (§'')"+i ^ x T be the linear interpo- 
lation operator given by 

Wn{u,v){t) := {[nt] + l~nt) (u[„t], U[„t]) -|- {nt - [nt]) (u[„t]+i, U[„t]+i), t £ [0, 1] 

where u — {uo^ui, ...,Un), v = {vo,vi, ...,Vn) and [z] denotes the integer part of z. 
Set 

(2.9) sup EF(W„(A/,(M))), 

we denote by E the expectation with respect to the underlying probability measure. 
The following theorem which will be proved in Section [5] is the main result of the 
paper. 

Theorem 2.2. For any e > there exists a constant Ce — Cc(j^) which depends 
only on the distribution v such that 

(2.10) |V;r - V| < Cen'"^/^ Vn e N. 

Furthermore, if the function F is bounded, then there exists a constant C — C(j') 
for which 

(2.11) \V,'^ ~V\<Cn-'^/^, VneN. 
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Next, we describe a dynamical programming algorithm for V.^ and for the op- 
timal control, which in general should not be unique. Fix n e N and define a 
sequence of functions J^'" : (M'^)'=+i x (S'*)'=+i -> R, fc = 0, 1, n by the backward 
recursion 

(2.12) J!;:^''^{uQ,ui,...,Un,VQ,vi,...,Vn)^F{W„{u,v)) and 

J^'"(uo,"i, ...,Ufe,wo,wi, ...,Ufe) = 

sup^e\/D/K<J -^fe+i ("o,ui, ...,Uk,Uk + ^,vo,vi, ...,Vk,Vk + v) ^vix) 
for fc = 0, 1. 
From (|2.ip and (|2.6p it follows that there exists a constant such that 

k 

J^'" (mo, .■.,Uk,vo, ...,Vk) < iJexp ((i72+l)^(||wi|l + lki|l)). Vfc,uo, ...,Ufc,t'o, ...,t;fc. 

Fix k. By applying (|2.6p again we conclude that for any compact sets Ki C M"* 
and K2 C S"^ , the family of random variables 

-'fc+i \ UQ,...,Uk,Uk + ^^7^,uo,...,Wfc,Wfe + IT J , 

7 e Vd, Wo, ■•■,"fc e i^i, vo,...,VkeK2 

is uniformly integrable. This together with the fact that the set D is compact gives 
(by backward induction) that for any k, the function J^'" is continuous. Thus we 
can introduce the functions /i^'" : (E'')'=+i x (S'')'=+i ^ a/D, fc = 0, 1, ...,n- 1 by 

(2.13) /i^'"(uo, ...j-UfejWo, •••,Wfe) = 
argmax^gyg/jj^ J^:^^ (^wo, Mi, Mfe, Ufc + ^,Vq,Vi, ...,Vk,Vk + ^) (ii/(a;). 

Finally, define by induction the stochastic processes {M^'^}^^^ and {^^'"l^^o' 
with values in M'' and S'^, respectively by A/q'" = 0, Nq'^ = and for fc < n 

(2.14) n;:-i\ = n;:^- + 1 (/ir«'", ^4''", K'" . 

and = Mr+^/»r(<''\--.^r-<"\--.<'")n+i. 

Observe that Af'" G and N"-^^ = (M'''"). From the dynamical programming 
principle it follows that 

(2.15) V;^ = Jo'"(0>0) = EF (W„(M"'", (iVr-"))) . 

In the following theorem (which will be proved in Section 5) we provide an explicit 
construction of e-optimal measures for the G-expectation which is defined in p.3p . 

Theorem 2.3. Let (fi^, J"^,P^) be a complete probability space together with 
a standard d- dimensional Brownian motion {Wt}te[Q,i] '^'^'^ natural filtration 
= (t{VI^(s)|s < t}. Consider the standard normal distribution Vg = A/'(0,/). 
For any n G N, let fn ■ (M**)" — > M &e a function which is satisfying fniYf, Yn) — 
Mn"'^ , where Yf,...,YJ^ are i.i.d. and C{Yf ) = i^g. Observe that fn can be calcu- 
lated from Ii2.12]) - ^2.14\ l- Define the continuous stochastic process {-/^"Ij^q by 
(2.16) 
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where denotes the expectation with respect to P^. Let Pn he the distribution 
of M" on the canonical space ^l. Then P„ G "Pd, and for any e > there exists a 
constant such that 

(2.17) V <EnF{B,{B))+C,n'-^/^, Vn 

where denotes the expectation with respect to P„ . // the function F is bounded 
then there exists a constant C for which 

(2.18) V < EnF{B,{B)) +Cn-^^^, Vn. 

3. The main tool 

In this section we derive a strong approximation theorem (Lemma 13. 2p which is 
the main tool in the proof of Theorems I2.2H2.3I This theorem is an extension of 
the main result in |11| . 

For any two distributions i'i,h'2 on the same measurable space {X,B) we define 
the distance in variation 

(3.1) p(l/i,i/2)= sup|i/i(B)-z.2(S)|. 

see 

First we state some results (without a proof) from [11] (Lemmas 4.5 and 7.2 in 
[11] ) that will be used in the proof of Lemma 

Lemma 3.1. 

i. There exists a distribution fi on R'' which is supported on the set (—1/2,1/2)'' 
and has the following property. There exists a constant Ci > such that for any 
distributions vi , V2 on which satisfy 

(3.2) /jjrf xdvi{x) ~ Jjga xdv2{x) and for I < i, j < d 

JiRd x'x^di^i{x) = x''x^di^2{x) 

we have 

(3.3) p{vi* ii,V2* fi) <Ci[ I \\x\\^dvi{x)+ I \\x\f'dv2{x)] 
where v * pL denotes the convolution of the measures v and /i. 

a. Let (O, P) be a probability space together with a d-dimensional random vector 
Y , a m -dimensional random vector Z (m is some natural number), and a random 
variable a which is distributed uniformly on the interval [0, 1] and independent of 
Y and Z . Let v be a distribution on and let v he a distribution on M™ x 
such that y{A x W^) = P[Z e A) for any A € B{W^), i.e. a marginal distribution 
of on M™ is equals to C{Z). There exists a measurable function $ — '^u,i>,c(z,Y) ■ 
X M'^ X [0, 1] ^ K'' X M'^ such that for the vector 

(3.4) {U,X):^^Z,Y,a) 

we have the following: C{U) — v, C{Z,X) — v, U is independent of X,Z and 

(3.5) PiU + X ^ Y\Z) = p{C{U) * C{X\Z),C{Y\Z)). 

Now we are ready to prove the main result of this section. For any stochastic 
process Z = {Zk}^^^ we denote AZ^ := Zk — Zk-i for fc > 0, where we set Z-i = 0. 
Fix n e N and consider a d-dimensional martingale {MkYl^Q with respect to its 
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natural filtration, which satisfies Mq = 0. For any 1 < fc < n, let 0fe : (R'')*^ — )• S"* 
be a measurable map such that 



(3.6) x/A(M)fc = (AA4AM;|a{Mo, Mi, Mfe„i}) = 

0fe(AMo,AMi,..., AMfe_i), 

where {(A/)fe}]J^g is the predictable variation (S"^ valued) of Af and the symbol 
denotes transposition. We assume that there exists a constant H for which 



(3.7) E(||AMfc||3|a{Mo,...,Mfc_i}) + ||VA(M)fe||3 < H, a.s. Vfc. 

Lemma 3.2. Let v a distribution on such that 

(3.8) /jjrf xdv{x) = 0, x'x^diy{x) = (5,j ViJ < d 

and Jjjj^ ||a;||^(ii^(x) < oo. 

For any Q > its possible to construct the martingale {Mk}^^^ on some probability 
space {Q,,T,P) together with a sequence of i.i.d. random vectors Yi,...,y„ with the 
following properties: 
I. C{Yi) = V. 

a. For any k, the random vectors Mi, ...,Mk-i are independent ofYk. 

Hi. There exists a constant C2 = C2(i^) which depends only on the distribution v 

such that 



CoHn 



(t/fc,Xfe) = *^,,,,p,(Xo,...,Xfc_i,t/fc_i + -<^fc(eXo,...,eXfc_i)Ffe,afe), l<k<n 



(3.9) P 1^ max^ | \M, ~ ^ ^^^M),Y, 1 1 > 9 j < ^3 . 

Proof. Fix > 0. For any k let Vk be the distribution of the random vector 
^(AA/q, AA/fe). Let {VL,P,P) be a probability space which contains a se- 
quence of i.i.d. random vectors Yi,...,y„ such that C{Yi) = a, sequence of 
i.i.d. random variables ai,...,a„ which are distributed uniformly on the interval 
[0, 1] and independent of Yi, Yn, and a random vector C/q which is independent 
of Yi, K„, ai, a„ and satisfies £(t/o) = where the distribution /i is defined 
in the first part of Lemma [3.11 Define the sequences {^i}"^o {^^4}"=! 
following recursive relations, Xq = and 
(3.10) 

1 

where Vk is the distribution of {Xq, Xk-i,Uk-i + ^4>k{QXf), ...,QXk-i)Yk). 
From the definition of the map ^E* it follows (by induction) that £(9Xo, QXn) — 
£{AAIo, AMn). We conclude that the the stochastic process QY^i^Q Xi, < 
< n is distributed as {Affej^^g, and so we set, 

k 

(3.11) Mk = e ^ Xj, 0<k<n. 

i=0 

Let 1 < fc < n. From p.lOI) - (|3.1ip and the fact that Yk is independent of 
Yi,...,Yfc_i, ai,...,afc_i it follows that Yk is independent of Mq, Mk-i. Thus 
in order to complete the proof, it remains to establish p.9p . Set 

(3.12) ^Uk+Xk- Uk-i - ^(pki^Xo, eXk-i)Yk, and pk{xo, ■■■,Xk-i) 
= P{Skj^O\Xo^xo,...,Xk-i^Xk-i), xo,...,Xk-i eR'^ l<k<n. 
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From the properties of the map ^ it fohows that for any k, Uk is independent of 
Xq, ...,Xk and C{Uk) = This together with p.Sp and p.lOp gives 

(3.13) Pk{xo, ...,Xk-i) = p{C{Xk\Xo = .To, ...,Xk-i = Xk^i) * 
Ci^<|){exo,...,exk-l)Yk)*^l) xo,...,Xfc_i e l<k<n. 

From (I33D, (ISS1)-(EI3, (IHUl) and (IXTa 

(3.14) pfc(xo,...,Xfc_i) < a;o,...,a;fc_i e R'^, 1 < /c < n 

for some constant C2 = ^2(1^) which depends only on the distribution ly. From 
p.lip - p.l2p . p.l4p and the fact that maxo<A:<„ \ \Uk\\ < | a.s. we obtain 

P (maxi<fc<„ llMfe - ^^ti ^/A{M)~Yj\\ > o) = 

P (maxi<fc<„ p4 - Ej^o <^j(AMo, AAf,)r,+i|| > o) = 

P (maxi<fc<„ e|| Eti S^ + Uo-Uk\\>e) < ELi 7^ 0) < ^ 
and we conclude the proof. □ 

4. Auxiliary lemmas 

In this section we derive several estimates which are essential for the proof of 
Theorem 12. 2H2.3I We start with the following general result. 

Lemma 4.1. Let {M^jj^Q be a one dimensional continuous martingale which sat- 
isfies '^^^'P' < "H a.s. for some constant H. Consider the discrete time martingale 
Nk = Mk/n, < k < n together with its predictable variation process {{N)k}^^Q 
which is given by {N)o = and 

k 

{N)k = ^E((A7V,)'k{^^o, ...,iV^-i}), l<k<n. 

1=1 

There exists constants 03,64 (which depend only on % such that) 

C3 



(4.1) E max max \Mt - Nk\ ] < 

\0<fe<n-l fe/ri<t<(fc+l)/n / 71 

and 

(4.2) e( max max \{M)t - {N)kn < 

\0<k<n-lk/n<t<{k+l)/n J 

Proof. From the Burkholder-Davis-Gundy inequality it follows that there exists a 
constant ci such that 

(4.3) E (maxo<fe<„_i riia-Xk/n<t<(k+i)/n \Mt - Nk\'^) < 



ELo ^ (niaxfc/„<i<(fc+i)/„ \Mt - Mfc/n]") 



< 



this completes the proof of (|4.ip . Next, we prove (|4.2p . Define the optional variation 
of the martingale {Nk}^^^ by [N]o = and 

k 

(4.4) [7V]fe = ^(AiV,)', l<k<n. 
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From the relation E{A[N]k\a{No, ...,Nk-^i}) = A{N)k and the Doob-Kolmogorov 
inequaUty we obtain 

(4.5) E (maxo<fe<„ \[N]k - {N)k\^) < 4E (|[iV]„ - (iV)„|2) ^ 

4E (I ELi A[7V], - A(iV),|2) = 4E:Li E (|A[iV], - A(7V),p) < 

4 ELi E((A[iV]02) = 4 E (|Af,/„ - M(,_i)/„|4) < 4c^ 

where the last inequality follows from the Burkholder-Davis-Gundy inequality. 
Next, observe that 



(4.6) [N]k =Nl-2 V N,{N,+^ ~ N,) = - 2 f ^ N[nt]dMt, 1 < fc < 

From the Doob-Kolmogorov inequality and Ito's Isometry we get 

(4.7) E (supo<„<i |/;(Mi - 7V[„,])dMtf ) < 

4E (l J^{Mt - iV[„t])dAft|2) = 4E (/o (Mt - N[,,t])'' d{M) t) < 



AHE (maxo<fe<„-i raa.Xk/n<t<{k+i)/n \Mt - iV^p) 



< 



the last inequality follows from (|4.1I) and Jensen's inequality. From (I4.6p - (|4.7|) and 
the equality 2 J^^'' MtdMt ^ ~ {M)k/n it follows that 

^l<fe<n ' J — 

This together with (|4.5p and the inequality (a + + c)^ < 4(a^ + 5^ + c^) yields 

(4.8) E (maxo<fe<„_i maxfc/„<t<(fe+i)/„ \{M)t - (iV)fe|2) < 

^ + 4E (maxi<fc<„ - (M)fc/„|2) + 4E (maxi<fc<„ | [TV]^ - (7V)fc|2) 

and the proof is completed. □ 

Next, we apply the above lemma in order to derive some estimates in our setup. 

Lemma 4.2. Let n G N and P G Vd- Consider the d-dimensional martingale 
Nk = -Bfc/n, < k < n together with its predictable variation {(-/V)fc})J^Q, under the 
measure P. There exists a constant C5 (which is independent of n and P) such that 

(4.9) Ep {\\Wn{N) - B\\^) < ^ 
and 

In the equations ^.9^ and ^.10\ l, Wn is the linear interpolation operator which is 
defined on the spaces (K<i)"+i and (8'')"+-'^, respectively. 

Proof. Inequality (|4.9p follows immediately from (|4.ip and the relation 

d 

||W„(A^)-i?||<2^max max \Nl-Bl\. 

^ l<k<n k/n<t<{k+l)/n 



(4.10) Ep{\\Wr.m)-{B)\\^)< 
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Next, we prove (|4.10p . For any 1 < i,j < d denote by {NY^-' and (B)]'-' , the i— th 
row and the j— th column of the matrices {N)k and {B)t, respectively. Notice that 
{B)l'' = l{{B'^ + B^)t-(B^)t~{B^)t) and (7V)^'^' = 1{{N^ +N^),- (N^),^ (N^),). 
Thus (I4.10p follows from (|4.2I) and the inequality 



|W„(W)-(i?>||<2^^ max max \{N)]f - {B)'^^^ 

— ' — ' 0<fc<n— 1 k/n<t<(k+l)/n 



. , . , - -k/n<t<{k+l)/ 

1=1 j = l 



We conclude this section with the following technical lemma. 

Lemma 4.3. Let A > 0. Then we have: 
i. 

(4.11) sup Epexp{A sup ||B4||) < oo. 

PeVo o<t<i 



□ 



ii. Let n gN and v he a distribution which satisfies \2..^) - \2^) . Consider a filtered 
probability space (f^, {-T^fcl^^o' ^) together with a sequence of i.i.d. random vec- 
tors Yi, Yn which satisfy C{Yi) = v. Assume that for any i, Yi is Ti measurable 
and independent of Ti-\. Let be a d~dimensional stochastic process of 

the following form: Mq — and 



(4.12) Mfe = -Vi^7,y„ l<k<n 



where for any i, •fi is J-i-i measurable random matrix, which takes values in \/D. 
There exists a constant Cg ( which may depend on A and v ) such that 



(4.13) exp A max ||Mfe|| < Cq. 

\ 0<k<n j 

Proof, i. Let P E Vri- From the Novikov condition it follows that for any 1 < i < d 
and a e M, exp (^aBl - = 1. Thus 

Ep (exp(a|Bl|)) < Sp(exp(aSi)) + Sp(exp(-aBl)) < 2 exp (yl|D||) 

where ||D|| = suppg^ | |. This together with the Cauchy-Schwartz inequality 
completes the proof of (|4.1ip . 

ii. Consider the compact set K :— {x <E R'^ : \\x\\ < ||a/D||}. Clearly, the rows 
of the matrices 7j, 1 < J < n are in K. Fix 1 < i < d and consider the i— th 
component of the process M, namely we consider the process {Mq, ...,M^). From 
()4.12p we get that for any a e M 

E fexp(a(Af^ - Ml_,))\A-i) < sup V, 



yeK 

where ipi, is the function which is defined below (|2.5I) . This together with (12. 6p 
gives 

(4.14) E (exp(aAf;)) < sup sup V" f ) < oo. 
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From the inequality Eexp(|aM,\|) < E (exp(aAf^))+E (cxp(— aAf^)) and the Cauchy- 
Schwartz inequality it follows that there exists a constant C2 (which may depend 
on A and v) such that 

(4.15) E(exp(yl||Af„||)) <C2. 

Finally, since for any i the process M^, fc < n is a martingale with respect to the 
filtration {J^k}k=Q we conclude that the stochastic process {exp(A||Mfc||/2)}5!^Q 
is a sub-martingale and so, from (|4.15p and the Doob-Kolomogorv inequality 
Eexp(Amaxo<fc<„ ||Affe||) < 4c2 and the proof is completed. □ 



5. Proof of the main results 



In this section we complete the proof of Theorems I2.2H2.3I Let v he a distri- 
bution which satisfies p.5p - (|2.6|) . Fix e > 0. We start with proving the following 
statements 

(5.1) >V~ C,n'-l/^ Vn e N 
and for a bounded F 

(5.2) V,"^ >V- CtT^I^, Vn e N. 

Choose n e N and S > 0. There exists a measure Q e for which 

(5.3) V <5 + EqF{B,{B)). 

Consider the stochastic process = B^j^, Q < k < n together with its predictable 
variation {{N) k}^^^. From the fact that D is a convex compact set we obtain that 
there exists a sequence of functions V'j '■ (K'')-' ~> "^/D, ^ < j <n such that 

(5.4) ^/^Wk = ( AiVfc ATV; I a{iVo , A^i , . . . , iVfe- 1 }) = 

^7^fc(7Vo,...,iVfc_i), Vfc a.s. 

From the Burkholder-Davis-Gundy inequality it follows that there exists a constant 
C3 for which 

(5.5) SQ(||AiVfe||3|a{Aro,...,iVfe_i}) <c3n-3/2, Vfc a.s. 

By applying p.ip . Lemmas I4.2H4.3I and Cauchy-Schwartz inequality we get 

(5.6) Eq\F{B,{B)) - F{yVn{N),yVn{{N)))\ < an-^l^ 

for some constant C4 (which depends only on the distribution v). From (j5.5p and 
Lemma 13.21 we obtain that there exists a probability space (f2,J-", P) which con- 
tains the martingale N, a sequence of i.i.d. random vectors li, ...,y„ and satisfies, 
C{Yi) = V, for any k the random vectors Ni^ Nk-i are independent of Yk, and 

fc 



for some constant C5 which depends only on the distribution ly. Denote Mk — 
Ej=i VA(iV)jY,-, 1 < fc < n and ^ = {maxi<fc<„ ||7Vfe - Mfe|| > n"!/®}. From 
(|2.5|) and the fact that Ni, Nk-i arc independent of Yk we obtain that M is 
a martingale, and (Af) = (A^). Thus from ((5?7|) . Lemma lOl the Markov 
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inequality and the Holder inequality (for p ~ ^nd q — we get that there 
exists constants cg, cy which depend on e and v such that 

(5.8) E\F{WniN),Wn{{N))) - ^^(W„(Af), W„((A/)))| < 

HiE[ cxp(i?2(niaxi<fe<„ ||M||fc + maxi<fc<„ ||iV||fc + 2||D||)) x 



{n~^/^+lAi\\Wnm\ + \\WniMm)j < 

< C6(n-i/8 + P(A)T^) < cyn-^-i/s 

where we set = 1 if an event A occurs and = if not, and E denotes the 
expectation with respect to P. If the function F is bounded, say F < R, then we 
have 

(5.9) ^|F(W„(iV),W„((iV)))-F(W„(M),W„((M)))| <i?P(A) + i/in-i/8 
x£;(exp(iJ2(maxi<fe<„ | |M| |fe + maxi<fe<„ \\N\\k + 2||D||))) < cgn^i/s 

for some constant cg which depends only on v. Since S > was arbitrary, then in 
view of ([5?^ and ((5^ - ((01) we conclude that in order to prove ([5TT 1) - ([0)) 

it remains to establish the following inequality 

(5.10) y,r > EF{WniM), W„((M))). 

Define a sequence of functions Lk : (R'')'=+i x (§'')'=+i ^ R, fc 0, 1, ...,n by the 
backward recursion 

(5.11) L„(uo,...,u„,wo,...,u„) =i^(yV„(M, w)) and 

Lk{uQ,...,Uk,VQ,...,Vk) = 

ELk+i{uo, + -i='0fc+i(uo, ...,Ufe)Yfc+i,wo, —,Vk, 

Vk + ^ipl+i{uQ,...,Uk)) for /c 0, 1, 1. 

From the fact that Yfc+i is independent of Yi, ...,Yk, Ni, Nk-i it follows (by 
backward induction) that for any fc, 

(5.12) E{F{Wn{M),Wnm)))W{Nu...,Nk-uYi,...,Yk}) = 

Lk (Mo,...,A4,(7V)o,...,(iV)fe). 

Finally, from (|2.12p . (I5.1ip - (|5.12p and the fact that tpk takes values in \/D for any 
k, we obtain (by backward induction) that Lk < ./fc'", k < n, and in particular 

(5.13) = ■/o'"(0'0) > Lo(0,0) = ^F(W„(M),W„((M))). 

This competes the proof of (|5.ip - (|5.2p . Next, fix n e N, a distribution i' which 
satisfies (I2.5I) - (|2.6I) and consider the optimal control Af^'" which is given by (|2.12l| - 
(12.141) . By applying Lemma [3.21 for the standard normal distribution i^g it follows 
that there exists a probability space {ft,J^, P) which contains the martingale M"^'", 
a sequence of i.i.d. standard Gaussian random vectors (d-dimensional) Yf ,...,Y^ 
such that for any k the random vectors M^"'", M^'"^ are independent of Y^ , and 

(5.14) P f maxJlM,^'" - ^ ^ A{M'^^-),Yf\\ > n-'A < cn^^/^ 
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for some constant cg. Denote Mk — J2'j=i \/A(M'^'")jY^^, 1 < k < n. Observe 
that (M) = (A/'''"). Thus by using similar argument to those as in (I5.8p - (|5.9p we 
obtain that there exists constants cio,cii such that 

(5.15) \EF{Wn{M),Wni{M))) ~ V^l < 
E\F{WniM),Wn{{M))) - F(W„(Af''>"), W„((M'^'")))| < cion^-i/^ 

and if the function F is bounded, 

(5.16) \EF{Wn{M),Wni{M))) ~ V;rl < 
E\F{WniM),Wni{M))) - F(W„(M'''"), W„((Af'''")))| < cnn-i/s. 

By applying similar arguments to those as in (j5.1ip - (|5.13l) we conclude that 

(5.17) O = ■/o""(0, 0) > EF(WniM),Wni{M)))- 

Next, let z/j : (M"*)*^ — > \/D, l<fc<n — Ibea sequence of functions such that 
for any 1 < fc < n - 1, zUYf , ...,Y,^) = /i:-"(AC-, M^-", <-"), 
where the terms M'''"", A^^s'" are given by (|2.12p - (l2.14p . From the martingale 
representation theorem if follows that the martingale Af„ which is defined by (|2.16p 
equals to 

Af," = C"(0,0)W^t+It>i/„x 

Il/n Zlnu]iVnWi/^, ^A^(W^2/n - W^l/n), %/«(W^[n«] " -1 ))rfW^« , t € [0, 1] 

and so we obtain that P„ e Vu- As in (15. 6p we have 

(5.18) i;„|F(B, (B)) - F{Wn{N),Wn{{N)))\ < c^n'^'^ 

where — B^/n, < k < n. Finally, observe that the distribution of N under P„ 
equals to the distribution of the martingale Af "". Thus from (|2.15p and (jS.lSp we 
conclude that 

V > Ep,^F{B, (B)) > O - C4n~i/^ 
This together with (|5.ip - (|5.2p and (I5.15l) - (|5.17p completes the proof of Theorems 

[2:2m:31 □ 
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